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Abstract 

The mean intensity of planetary nebulae with an expanding atmosphere is modeled 
by considering dusty and dust-free atmospheres. The bulk matter density is deter- 
mined from the adopted velocity field through the equation of continuity. The gas is 
assumed to consist of hydrogen and helium and the gas-to-dust mass ratio is taken 
to be 3 x 10~ 4 . The Rayleigh phase function is employed for atomic scattering while 
the full Mie theory of scattering is incorporated for determining the dust scattering 
and absorption cross-section as well as the phase function for the angular distri- 
bution of photons after scattering. It is shown that in a dust free atmosphere, the 
mean intensity increases with the increase in the expansion velocity that makes the 
medium diluted. The mean intensity profile changes significantly when dust scat- 
tering is incorporated. The increase in forward scattering of photons by the dust 
particles yields into an increase in the mean intensity as compared to that without 
dust. The mean intensity increases as the particle size is increased. Thus it is shown 
that both the expansion of the medium and the presence of dust play important 
role in determining the mean intensity of a planetary nebulae. 
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1 Introduction 



There are several models of planetary nebulae which describe that the ioniza- 
tion structure and other physical characteristics (Hjellming 1966, Rubin 1968, 
Harrington 1968, 1969, Kirkpatrick 1970, 1972, Koppen 1979, 1980 ). They 
calculated the ionization structure by solving simultaneously, the equations 
of radiative transfer and ionization equilibrium for several elements such as 
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H, He, C, N, O, Ne, Mg, S, Ar etc. in different stages of ionization together 
with the energy equilibrium equation. The ionization equilibrium equation re- 
quires the use of the radiation in the form of mean intensity. Normally mean 
intensity obtained by the 'on the spot approximation' is utilized in evaluating 
the formation of hydrogen Ly a line in dusty, expanding planetary nebulae. 
Infrared observations of planetary nebulae have established in the presence 
of dust (like graphite, amorphous carbon, silicate, and iron) in these objects 
and a comprehensive review work up to 1982 has been presented by Barlow 
(1983). Pottasch et al (1984), Iyengar (1986), Zijlstra et al (1989), Ratag et 
al (1990), Amnuel (1994). Pottasch et al. (1984) found infrared excesses in 
high dust temperature nebulae. They concluded the possibility of dust be- 
ing mostly heated by the radiation of the central star on the long wave side 
since the nebulae are young. Ratag et al. (1990) suggested that the infrared 
excesses was due to the dust contents or heating by the interstellar radia- 
tion field. Amnuel (1994) favored the latter explanation. Peraiah and Wehrse 
(1978) found that the mean intensities of the radiation field are far different 
from those obtained from the 'on the spot approximation'. The expansion in 
the nebulae would certainly modify the radiation field and this radiation field 
in tern changes the ionization structure in the nebulae. However, a complete 
and rigorous analysis by incorporating relevant theory for dust scattering in 
an expanding medium is not explored. In this paper we investigate the changes 
that are produced on the radiation field by the radially expanding gases and 
dust. For this purpose, we develop the solution of radiative transfer equation 
in the expanding nebulae. We assume both gas and dust with hydrogen and 
helium as the components of the gas. In section 2 we present the relevant 
radiative transfer equations that describes the radiation field under consider- 
ation. In section 3 we present the absorption cross section for atoms and dust. 
The adopted velocity law, the density distribution and the dust parameters 
are presented in section 4. The numerical method is described in section 5. 
The results are discussed in section 6 followed by our conclusion in section 7. 



2 Radiative Transfer Equations 



In a spherically symmetric, expanding, dusty planetary nebulae, the equa- 
tions of radiative transfer can be written (Peraiah and Grant, 1973; henceafter 
PG73) as : 
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for the inward- going rays. 

The quantity [/ (r, /x) is given by, 

[/ (r, /x) = 4nr 2 I(r, /x) 



(3) 



where J(r, /x) is the specific intensity of the ray making an angle cos -1 /x with 
the radius vector at radial points r. We have restricted /x = cos# to lie in the 
interval (0,1). Similarly, 



B g , d (r) = 4itr 2 b g4 (r) 



(4) 



where the quantities b 9td (r) are Planck functions for gas and dust respectively. 
ujg and oj d are the albedo's for single scattering for gas and dust. The quantities 
o g and a d are the absorption cross-sections for gas and dust. 



3 Absorption Cross-section for Gas and Dust 



For the absorption cross-section of hydrogen and helium we have (Hummer 
and Seaton 1964) 
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for helium. Here z/ is taken to be 3.289 x 10 15 Hz which is hydrogen Lyman 
limit. 

The calculations are being done in the rest frame and therefore, we have to 
modify the absorption cross section an and a He in the rest frame by taking into 
account of the Doppler effect. The equation (5) and the (6) can be rewritten 
as (Landau and Lifshitz 1975) 
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velocity of light. We have restricted our calculations up to a maximum fre- 
quency equal to v = 4 x u . From equations (7), (8), (10) and (11) it is clear 
that 



a g (r, /i, v) = a+(r, v) = a%(r, v) + a% e (r, v) 
a g (r, -//, u) = <j-(r, v) = a H {r; u) + ajj e {r; v) (13) 

We have employed Rayleigh phase function that governs the angular distribu- 
tion of photons after scattering with the atoms and is given by (Chandrasekhar 
1960) 



P(r, w') = I [3 - + (2 M 2 - l)// 2 ] (14) 



For dust absorption and scattering we have employed the Mie theory of scat- 
tering (van de Hulst 1957). The extinction and scattering cross-section for 
dust scattering are given by 
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where a is the radius of dust particles and x = 2ira/\. The Mie coefficient 
a n and b n are in general complex function of the refractive index of the dust 
and of x. We have calculated the dust extinction and scattering cross section 
and hence the albedo Ud numerically for specific values of the dust size and 
refractive index. We have employed the full Mie phase function for the angular 
distribution of photons after scattering with dust. 



4 Model Parameters 



We have adopted the usual power law for velocity profile which is given as 
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V(r)=V + V t (l-^)\ (17) 

where V is the initial velocity at the inner radius R in = 9.2568 x 10 16 cm, and 
V t is the terminal velocity at R out = 2.221632 x 10 17 cm taken from Osterbrock 
(1974, page 135) . We fix the initial velocity Vq at 0.1 kms^ 1 . The terminal 
velocity V t is considered to be 10, 20, 30, 50, 100 and 300 kms -1 . 

The distribution of total density comes from the equation of continuity and is 
written as 



"< r > = i^W) < 18 » 

where M is the rate of mass loss. Since the particle number density in planetary 
nebulae ranges from 0.1 to 100 cm~ 3 , we have set M = 10~ 6 M yr _1 . The 
number density of He is taken to be 0.1 times that of H. 

The dust-to-gas mass ratio for IC 418 and NGC 7662 are derived by Hoare 
(1990) to be 6 x 10~ 4 and 3 x 10~ 4 respectively. In our model we have taken 
dust-to-gas mass ratio as 3x 10~ 4 . It is generally believed that dust in planetary 
nebula is carbon based one such as graphite. However, observation at infrared 
reveals the presence of silicates (Stasinska & Szczerba 1999; Pottasch 1987). 
Thus, we consider silicate type grains to demonstrate the effect of forward 
scattering. Now both the real and the imaginary refractive index for silicate 
grains are function of frequency and they depend on the grain composition. 
The real part of the refractive index for magnesium silicates such as Forsterite 
and Enstetite varies from 1.55 to 0.8 while the imaginary part varies from 0.85 
to 0.3 (Scott & Duley 1996) in the frequency interval we have considered in 
the present work. For the sake of simplicity, in the present work we have taken 
an average value of the refractive indices for all the freequency points. The 
real part of the refractive index is taken to be 1.2 while the imaginary part 
is set at 0.5 for all frequency points. We have assumed same kind and size of 
spherical dust particles with radius 0.1 and 0.5 micron. It is worth mentioning 
that the Mie phase function reduces to that of Rayleigh when the ratio of 
the particle size to the wavelength is very small. Mie scattering is asymmetric 
in the sense that the forward scattering and the backward scattering are not 
the same while Rayleigh scattering is symmetric. The asymmetry parameter 
g =< cos 9 > is given by (Bohren & Huffman 1983) 
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where Q sca is the scattering efficiency and * denotes the complex conjugate. 
In fig 1, we present the asymmetry parameter for different frequency points 
with grain radii 0.1 jim and 0.5/zm. For Rayleigh scattering, g = because 
of symmetry. If the particle scatters more light toward the forward direction 
then g is positive while g is negative if the scattering is directed more to- 
ward the back direction. Although, the dust sizes considered is very small, 
the wavelength range that we are concerned is much shorter than the particle 
radius. As a consequence, g is positive in all frequency points considered in 
the present work. The numerical value of g very near to +1 implies that the 
particles scatter light mostly in the forward direction. For a given frequnecy 
point the asymmetry increases as the particle size increases. On the other 
hand, for a given particle size, the asymmetry increases as the frequency in- 
creases. However, the asymmetry parameter depends on the refractive index 
as well. With the increase in the refractive index, the absorbity of the particle 
increases. Fig 1 shows that g is almost constant with the particle radius 0.5 
fim. When the imaginary part of the refractive index is decreased from 0.5 
to as small as 0.02, g increases with the increase in frequency and saturates 
at much higher value. The asymmetry parameter for different size of graphite 
grains is presented by Draine (1985). The Mie scattering and extinction cross- 
sections for astronomical silicates have been discussed in details by Draine & 
Lee (1984). At large value of x, the extinction efficiency Q ext approaches to 
2.0. It is worth mentioning here that the dust opacity plays dominant role over 
the phase function in determining the emergent flux. 



5 Numerical Method 

The solution of radiative transfer in spherical symmetry (1) and (2) is devel- 
oped by using discrete space theory of radiative transfer (PG73). In general 
the following steps are followed for obtaining the solution. 

(i) We divide the medium into a number of "cells" whose thickness is less 
than or equal to the critical (r crit ). The critical thickness is determined on the 
basis of physical characteristics of the medium. r crit ensures the stability and 
uniqueness of the solution. 

(ii) Now the integration of the transfer equation is performed on the "cell" 
which is two-dimensional radius - angle grid bounded by [r n , r n+1 ] x i , i] 

where A^+i= X/jUi C k , j — 1, 2 . . . , J, where C k are the weights of Gauss Leg- 
endre formula. 

(iii) By using the interaction principle described in PG73, we obtain the re- 
flection and transmission operators over the "cell". 
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Fig. 1. The asymmetry parameter g with respect to the frequency. Top solid line 
represents the value of g for the particle radius 0.5 fim while bottom solid line 
represents that for particle radius 0.1 jim. The dashed line represents the value of 
g for particle radius 0.5 fim but with imaginary refractive index 0.02 

(iv) Finally we combine all the cells by the star algorithm described in PG73 
and obtain the radiation field. 



Now equation (1) and (2) can be rewritten upon \i integration 
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o; g and are the albedo's for single scattering. Here, C, M, a~ are the 
diagonal matrices of quadrature weights, angles and absorption cross-sections 
corresponding to the set of /i's. A + and A~ are the curvature matrices defined 
in PG73, and 
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where < fij < fij < 1. We discretize equations (19) and (20) by integrating 
these equations from r n to r n+1 . This gives us 
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where the quantity with subscript n + \ represents the average over the 'cell' 
bounded by the radii r n and r n+ \ (here, the 'cell' is defined as a shell with radial 
boundaries r n and r n+ \ and whose optical thickness is less than, or equal to 
the critical optical depth for which a stable and unique solution exists). Here, 



Ar n+ i = r n+1 - r n 



t + , =(7 + ,Ar n+ i (26) 
and Td is parametrized. The quantity p c is the curvature factor given by, 



Ar 



(27) 



where r n+ i = \{j- n + r n+ i) is suitable mean radius and / is the identity matrix. 
The quantities U^ +L are approximated by those at the boundaries r n and r n+ \ 
by (see PG73) 



(28) 



Substitution of equation (27) into equation (23) and (24) would give us the 
system of equations in the form of interaction principle and the result is written 

as 
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Comparing this with the interaction principle, 
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where t's and r's are the transmission and reflection operators. The quantities 
1 are the source vectors. By comparing equations (28) and (30) one obtains 
the reflection and transmission operators. These are given in the appendix. 



6 Results and Discussion 



As the purpose of this paper is to show the qualitative changes introduced by 
curvature and expansion. We have chosen Harrington's (1969) model of NGC 
7662 given in Osterbrock (1974, page 135). Few discrete points in the frequency 
grid are considered starting 3.289 x 10 15 Hz to about four times this frequency. 
We have divided this frequency interval into 20 points equally distributed along 
the grid. The mean intensity is calculated by using the relation 



1 ? 

Jv = 2 J I ^ dfM - 



(32) 



The results are presented in figure 2-5. Figure 2 shows the mean intensity with- 
out the effect of dust scattering and absorption. The mean intensity increases 
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Fig. 2. Mean intensity without dust particles for different frequency points. From 
bottom to top the curves represent the models with terminal velocity Vt = 10, 20, 
30, 50, 100 and 300 kms" 1 respectively. 

with the increase in the velocity. The change in the mean intensity is maxi- 
mum when the terminal velocity Vt is increased from 10 kms^ 1 to 20 kms -1 . 
Expansion of the medium yields the atmosphere diluted by decreasing the den- 
sity. Hence more photons can escape from the bottom of the atmosphere. As 
a result, the mean intensity increases with the increase in velocity. For a dust- 
free atmosphere, the opacity is determined by the absorption cross section of 
Hydrogen and Helium. The absorption cross section for Hydrogen varies with 
the cube of the frequency while that for Helium changes exponentially. For 
small expansion velocity the mean intensity falls slowly with the increase in 
frequency up to 10 16 Hz but then falls rapidly. This means the optical depth 
with small velocity is mainly governed by the change in the Hydrogen opac- 
ity. However, the Helium opacity starts dominating as the frequency increases 
from 10 16 Hz. When the velocity increases, the mean intensity increases with 
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Fig. 3. Mean intensity with dust particles of different sizes and for Vj=10 kms . 
Top solid line represent the model with dust particles of radius 0.5 fj,m while the 
bottom solid line represent the model with dust particles of radius 0.1 fim. The 
dashed line represents the model without dust but with the same terminal velocity. 

the increase in frequency and peaks at about 5.5 x 10 15 Hz and then falls 
rapidly. The change in the qualitative feature of the mean intensity with the 
increase in velocity is dictated by the change in the absorption cross-section 
while the increase in the intensity is due to the dilution of the medium. 

With the incorporation of dust grains, there are significant qualitative as well 
as quantitative changes in the mean intensities as can be seen from figure 3 
to figure 5. This is due to the fact that in a dusty atmosphere the optical 
depth is determined mainly by the dust. The absorption and scattering by 
dust are determined by the size and the refractive index of the grains. In 
the present investigation we have considered dust grain with sizes 0.1 jim 
and 0.5 \xm. Given the high temperature of the medium, grains with larger 
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Fig. 4. Same as figure 3 but with 14=100 kms -1 . 

size would not survive. Since the wavelengths under consideration is shorter 
than 0.1 fim, scattering by dust plays a significant role in determining the 
mean intensity. It is well known that Rayleigh scattering is symmetric so 
that the amount of forward and backward scattering is the same. But, Mie 
scattering is asymmetric and it makes the forward scattering more than the 
backward scattering as discussed in section 4. This asymmetry increases as 
the ratio between the wavelength and the particle size increases. In a dust 
free clear atmosphere, a significant amount of radiation gets back scattered 
contributing to the reduction in mean intensity. But, dust scattering makes a 
large number of photons scattered in the forward direction yielding an increase 
in the mean intensity as compared to that without dust. With the adopted 
values for the wavelength and the particle size, dust scattering dominates over 
dust absorption. As a result, increase in mean intensity is obtained when dust 
is incorporated in the atmospheric models. The mean intensity increases with 
the increase in praticle size as the asymmetry increases with the increase in 
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Fig. 5. Same as figure 3 but with 14=300 kms . 

particle size. It is worth noting that the mean intensity for two different sizes 
of dust particles is the same for smaller frequency or longer wavelength and 
the changes in the mean intensity increases as one moves from smaller to larger 
frequencies. 



7 Conclusion 



The mean intensity from a planetary nebulae with expanding atmosphere is 
modeled by solving the radiative transfer equations for multiple scattering. 
Both dusty and dust-free medium is considered. For the dusty medium, Mie 
theory of scattering by spherical dust grains is employed. In a dust-free at- 
mosphere, the mean intensity increases with the increase in the expansion 
velocity as the expansion of the medium removes matter causing a decrease in 
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optical depth. The entire intensity profile changes when dust scattering is in- 
corporated. Due to the increase in forward scattering of photon by dust grains, 
the mean intensity increases significantly in a dusty medium as compared to 
that in a dust-free atmosphere. Also, the mean intensity increases with the 
increase in particle size. Thus it is shown that both the expansion and the 
presence of dust in the atmosphere of planetary nebulae play important role 
in determining the mean intensity of the object. The important implication 
of the present work is that even the presence of very small dust grain would 
affect the spectrum of planetary nebulae at the ultra-violet and near optical 
region. Therefore observation at shorter wavelengths would provide important 
information on the properties of dust grains as well as the velocity field. 
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A Appendix 

The transmission and reflection operators as are given by, 

t(n + l,n)= t + [A + S ++ + r + - r - + ] 
f (n, n + l)=t~ [A~S + r - + r + -] 
r(n + 1, n) = 2t~r~ + A + M 

r(n,n+ 1) = 2t+r + ~A~M (A.l) 
and the source vectors are 



17 



£+ = (1 - ujg)T g t + [A + B + + r + -A~B-} 
IT = (1 - w fl )r fl r[A _ S- + r- + A + B + ] 



where 
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